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Abstract 

We use the AdS/CFT correspondence in a regime in which the field theory reduces to fluid 
dynamics to construct an infinite class of new black objects in Scherk-Schwarz compactified 
AdSti+2 space. Our configurations are dual to black objects that generalize black rings and 
have horizon topology 5"*~" x T" for n < Locally our fluid configurations are plasma 

sheets that curve around into tori whose radii are large compared to the thickness of the 
sheets (the ratio of these radii constitutes a small parameter that permits the perturbative 
construction of these configurations). These toroidal configurations are stabilized by angular 
momentum. We study solutions whose dual horizon topologies are X S\ S"* X and 

X T^ in detail; in particular we investigate the thermodynamic properties of these objects. 
We also present a formal general construction of the most general stationary configuration of 
fluids with boundaries that solve the d dimensional relativistic Navier-Stokes equation. 
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A Extrinsic curvature 



B Energy, angular momentum, and entropy of the ring and torus in 4+1 dimen- 
sions 

C Notation 



1 Introduction 

In this paper we attempt to understand horizon topologies and thermodynamics of blaek objects 
in arbitrary high dimensional Scerk-Schwarz compactified AdS spaces (SSAdS). The spectrum of 
black objects in more than 4 dimensions is extremely rich and consequently has drawn considerable 
interest recently [0, |^. As the construction of these exotic horizon topologies directly in gravity 
turns out to be technically difficult we study them in a somewhat indirect manner using the 
AdS/CFT correspondence |-|. 

This correspondence relates a theory of gravity in AdS(j+2 to a CFT in d + 1 dimensional 
spacetime. We consider the field theory obtained by Scherk-Schwarz compactification of this 
dual CFT, which consequently lives in d dimensions. This field theory has a first order confine- 
ment /deconflnement phase transition. This corresponds to a Hawking-Page-like phase transition 
in the bulk, for which the low temperature phase is the AdS-soliton and the high temperature 
phase is a large AdS black brane 

In the long wavelength limit, this field theory admits a fluid description where the dynamics is 
governed by the d dimensional relativistic Navier-Stokes equation. The effect of the Scerk-Schwarz 
compactification is only to introduce a constant additive piece to the free energy of the deconfined 
fluid 0]. Due to this shift, the pressure can go to zero at flnite energy densities, allowing the 
existence of arbitrarily large finite lumps of deconfined fiuid separated from the confined phase by 
a surface - the plasmaballs of [Q. Now by the AdS/CFT correspondence finite energy localized 
non-dissipative configurations of the plasma fiuid in the deconfined phase is dual to stationary 
black objects in the bulk. Thus, by studying fluid conflgurations that solve the d dimensional 
relativistic Navier-Stokes equation we can infer facts about the black objects in SSAdS(i+2 ^• 

Two important feature of the dual black object that one can infer from the fluid configurations 
are the horizon topology and the thermodynamics. The thermodynamics of the black object can 
be studied by simply computing the thermodynamic properties of the fluid conflguration - one 
integrates the energy density, entropy density etc. to compute the total energy, entropy etc. and 
the rest follows. 

The horizon topology can be inferred as follows. Far outside the region corresponding to the 
plasma, the bulk should look like the AdS-soliton. In this conflguration the Scherk-Schwarz circle 
contracts as one moves away from the boundary, eventually reaching zero size and capping off 
spacetime smoothly. Deep inside the region corresponding to the plasma, the bulk should look 
like the black brane. In this configuration the Scherk-Schwarz circles does not contract, it still 
has non-zero size when one reaches the horizon. It follows that as one moves along the horizon, 
the Scherk-Schwarz circle must contract as one approaches the edge of the region corresponding 
to the plasma. The horizon topology is found by looking at the flbration of a circle over a region 
the same shape as the plasma configuration, contracting the circle at the edges |7[ ^. We have 
provided a schematic drawing of this in fig.0. 

In the fiuid description the degrees of freedom includes the velocity field, and the 

temperature field, T{x), (we consider uncharged fluids dual to uncharged black objects; otherwise 
the degrees of freedom would also include the chemical potentials for those charges) . Now as we 
seek time-independent solutions, Lorentz symmetry allows us to consider fluid velocities of the 
form 

= j{dt+UjJa), 

where la are the killing vectors along the Cartan direction of the spatial rotation group, 7 is the 
normalization and the Wa are some constants. This along with the fact that our solutions are 
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AdS soliton Black brane AdS soliton 



Figure 1: Schematic description of the bulk dual of a plasmaball with some circle fibres indicated. 

non-dissipative forces the temperature field to be of the form 

T = iT, 

where T is a constant. With a simple thermodynamic argument we show that T is the overall 
thermodynamic temperature of the fluid configuration and are the thermodynamic angular ve- 
locities. Further we demonstrate that the equations of motion for non-dissipative time- independent 
solutions at the surface of the fluid configuration reduce to the condition 

I surface — fJ0, 

where a is the surface tension and is the trace of the extrinsic curvature of the fluid surface 
under consideration. The pressure V is related to the temperature T by the equation of state, 
so this provides a differential equation for the position of the surface. These conflgurations are 
parameterized by the temperature T and the angular velocities Wq. 

We then proceed to construct a class of fluid configurations whose surface is a solution of the 
above equation in a certain limit. In d spacetime dimensions the topologies of these configurations 
are 

^ ^ ^ 

n times 

where n satisfies 

ri = 0, for d = 3, 
d- 1 

n < — - — , for odd d greater than 3, 

n < , tor even d. 

- 2 ' 

These solutions are rotating in the plane in which the S^s lie, and for simplicity we turn off 
angular momentum along any other directions. In these configurations pressure in the radial 
direction of the ball is balanced by the surface tension. While along the radial direction of the 
S^s the centrifugal force balances the pressure (therefore rotation is essential in the plane of the 
S^s). We refer to the limit in which the (average) radius of the ball is small compared to the 
(average) radius along the S^s Q as 'the generalized thin ring' limit. The ratio of these two radii 
serves as the small parameter in the problem. To leading order in this parameter we find that the 
fluid configurations are exactly _b('*~i~") x T" (in contrast to merely having the same topology). 
The force balance conditions then relate the intrinsic fluid parameters (the temperature and the 
angular velocities) to the parameters of the fluid conflguration (the radius of the ball and the radii 
of the various S'^s). These fluid conflgurations are dual to black objects with horizon topologies 

^When there is more than one this radius refers to the magnitude of the vector which is obtained by the 
vector sum of the radii of the various S^s 
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g{d-n) rpn hence this provides an indirect proof of existence of such exotic horizon topologies 
of black objects in SSAdSd+2- This approach is reminiscent of (and inspired by) the black- fold 
approach of Q . 

Further it is possible to construct a well controlled perturbation theory about these generalized 
thin ring solution. This we demonstrate by explicitly computing the leading order corrections to 
the thin ring solutions in specific examples, namely the ring in 4 dimensions and the ring and 
the 'torus' (the one with the topology x T^) in 5 dimensions. We find that the leading order 
correction only appear at the second order in the expansion parameter (the small parameter 
described above). Also in these cases we explicitly compute the thermodynamic quantities (which 
are again correct up to second order in the expansion parameter) with which we construct the 
phase diagrams of these solutions within appropriate validity regimes. 

The rest of the paper is organized as follows. In §H, we review the general formalism we will 
use - the thermodynamics of the type of confining theory we consider in this paper, relativistic 
fluid mechanics with surface tension and the construction and thermodynamics of equilibrium 
configurations. 

In §1^, we discuss the construction of configurations in arbitrary dimensions in the 'thin ring' 
limit to first order in the expansion parameter. We do not take the expansion any further, as the 
coordinate system we use is ill suited to perturbation theory. 

We then switch to adapted coordinate systems and study three examples in detail. In §^ we 
study the ring solutions in four and five dimensions (topologies x and B^ x respectively). 
We construct the solutions and study their thermodynamics to second order in the expansion 
parameter. In we do the same for the 'torus' solution in five dimensions (topology x T^). 

In Appendix we compute some formulae for the extrinsic curvature of the surfaces. In 
Appendix we present expressions for the energy, entropy and angular momentum of the ring 
and 'torus' in five dimensions that were too large to include in the main text. We summarise the 
notation used in this paper in Appendix O. 



Hydrodynamics and overall thermodynamics of plasma- 
lumps 



In this section we review the general formalism we will use in this paper. In §2.1 we discuss 
the thermodynamic properties of the fluids we consider here, in S|2.2| we review relativistic fluid 



mechanics, in §2.3 we review the relativistic treatment of surface tension and in §E.4 we describe 



the general construction of equilibrium configurations. 



2.1 Thermodynamics 

A fluid with all conserved charges and chemical potentials set to zero satisfies 

p + V = sr, 

dp = Tds, (1) 
dV = s dT, 

where p, V, s and T are the local density, pressure, entropy density and temperature as measured 
in the rest frame of the fluid. Note that all intensive thermodynamic quantities can be written as 
functions of one variable which we will usually choose to be the temperature. Once we are given 
the pressure as a function of temperature, we can use ([T]) to determine the other quantities. 

The fluids that we will consider here - those obtained by the compactiflcation of a conformal 
theory with a gravity dual on a Scherk-Schwarz circle - have an equation of state of the form 

V = ^^{T''+' (2) 

Note that this notation is slightly different from that used in ^. The quantity a here differs 
from the one used previously by a factor of %■ We also have po — oTc'^- 
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When the parent conformal theory is A/" = 4 super Yang-Mills, we have a = ^ . 

As the confined phase has a free energy ~ 0{N'^), to leading order at large N we can treat it 
as the vacuum. The confining/deconfining phase transition occurs when the two phases have the 
same free energy. In our case, this is approximately aX T — Tc- At this temperature, the density 
is given by pc ^ {d + l)aTc'^- Note that po is not the critical density. 

2.2 Fluid mechanics 

Provided all length scales are large compared to the thermalisation scale of the fluid (which we call 
^mfp), each patch of the fluid is well described by equilibrium thermodynamics in its rest frame. 
The fluid is characterised by the velocity of these patches — described by a vector u'^ = ^{l,v) 
— and the intensive thermodynamic quantities in their rest frames — which can all be computed 
from the proper temperature T using the equation of state and the first law of thermodynamics, 
as in 



2.1 



The equations of fluid dynamics are simply a statement of the conservation of the stress tensor 

V^T'"' = 0. (3) 

This provides d equations for the evolution of for the d quantities H and T once we have expressed 
the stress tensor as a function of these quantities. 

2.2.1 Perfect fluid stress tensor 

The dynamics of a fluid is completely specified once the stress tensor and charge currents are given 
as functions of T and m^. As we have explained in the introduction, fluid mechanics is an effective 
description at long distances (i.e, it is valid only when the fluid variables vary on distance scales 
that are large compared to the mean free path Zmfp)- As a consequence it is natural to expand the 
stress tensor in powers of derivatives. In this subsection we briefly review the leading (i.e. zeroth) 
order terms in this expansion. 

It is convenient to define a projection tensor 

P^'" = g''" + uf'ur (4) 

p^!^ projects vectors onto the (d— 1) dimensional submanifold orthogonal to u^. In other words, 
P^"^ may be thought of as a projector onto spatial coordinates in the rest frame of the fluid. In a 
similar fashion, —u^u'^ projects vectors onto the time direction in the fluid rest frame. 

To zeroth order in the derivative expansion, Lorentz invariance and the correct static limit 
uniquely determine the stress tensor, charge and the entropy currents in terms of the thermody- 
namic variables. We have 

^perfect — pu- u- -r I r- , 

(<-^5 )perfcct — SU^ ^ 

where all thermodynamic quantities are measured in the local rest frame of the fluid, so that 
they are Lorentz scalars. It is not difficult to verify that in this zero-derivative (or perfect fluid) 
approximation, the entropy current is conserved. Entropy production (associated with dissipation) 
occurs only at the first subleading order in the derivative expansion, as we will see in the next 
subsection. 

2.2.2 Dissipation 

Now, we proceed to examine the first subleading order in the derivative expansion. In the first 
subleading order, Lorentz invariance and the physical requirement that entropy be non-decreasing 
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determine the form of the stress tensor and the current to be (see, for example, §§14.1 of pi[) 

dissipativc 

" r 



V S ydissipative 



where 

are the acceleration, expansion, shear tensor and heat flux respectively. 

These equations define a set of new fluid dynamical parameters in addition to those of the pre- 
vious subsection: C is the bulk viscosity, r/ is the shear viscosity and k is the thermal conductivity. 
Fourier's law of heat conduction q = — kVT has been relativistically modified to 

g^ = -KP^''(a,r + a.r), (s) 

with an extra term that is related to the redshifting of the temperature. 

At this order in the derivative expansion, the entropy current is no longer conserved; however, 
it may be checked 111 that 



TV^J^s ^'L^+Cd^ + 277a,,,a^^ (9) 



As and a^^^ are all spacelike vectors and tensors, the RHS of (|9|) is positive provided ry, C, k are 
positive parameters, a condition we further assume. This establishes that (even locally) entropy 
can only be produced but never destroyed. In equilibrium, J^' must vanish. It follows that, q^ , 
d and cr'"^ each individually vanish in equilibrium. 

For fluids with gravity duals, the shear viscosity is given by 77 = ^ ||l2|. We can estimate the 
thermalisation length of the fluid by comparing coefficients at different orders in the derivative 
expansion 

Up ^ - = 7^- (10) 

This length scale may plausibly be identified with the thermalisation length scale of the fluid. 
This may be demonstrated within the kinetic theory, where Zmtp is simply the mean free path 
of colliding molecules, but is expected to apply to more generally to any fluid with short range 
interactions. 

With the equation of state (|^), this is given by 

^-^'"^ nidT^+^ + T/^'Y ^^^^ 
As we will be restricting attention to temperatures close to Tc, we have ^mfp ^ 1/7^- 



2.3 Surfaces 

The plasma ball configurations we consider have a domain wall separating a bubble of the decon- 
fined phase from the confined phase. As the density, pressure, etc. of the deconfined phase are a 
factor of N'^ larger than the confined phase, we can treat the confined phase as the vacuum and 
the domain wall as a surface bounding the deconfined fiuid. 

At surfaces, the density of the fiuid changes too rapidly to be described by fluid mechanics. 
However, provided that we look at length scales much larger than the thickness of the surface, we 
can replace this region by a delta function localised piece of the stress tensor. 
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At these length scales, this stress tensor will depend on the direction of the surface, with 
dependence on its curvature being suppressed. 

In general, introducing a surface energy density ctb, a surface entropy density as and a surface 
tension cr, considerations similar to those leading to lead to 

= <T + Tas, 
da = —as dT. 

However, the surface tension was only computed &i T — Tc in [Q, so we will have to ignore its 
temperature dependence. As we can see above, this is equivalent to setting as = and aE = a. 

Let's describe the location of the surface by a function f{x) that is positive inside the fluid 
and has a first order zero on the surface: 

T^''^0{f)TZ, + 6{f)TZ,.^,- (12) 

At large length scales, as mentioned above, T^^^^^^^^ will only depend on the first derivative of / 
and no higher derivatives. 

If we demand invariance under reparameterisations of the function f{x) — > g{x)f{x), where 
17(2;) > 0, and that the surface moves at the velocity of the fluid 

«^5m/I/=o = 0, (13) 
the surface stress tensor is (see §2.3 of [||) 



TsZf.co = Vdfdf [aEU^'W- - a{g^^ - n^n^ + u^u^)\ , (14) 

where = —df_if /y/df -df is the normal to the surface. Note that (9^/)Tsurfacc = 0. If we take 
the surface tension to be constant, as above, we get 



TsZ,.co = -<^h^''Vdrdf, (15) 

where Ti^^ = g^i, — n^rii, is the induced metric of the surface. The factor of y/df-df also has 
a simple interpretation: suppose we use a coordinate system where / is one of the coordinates. 
Then 

which provides the correct change of integration measure for localisation to the surface. If we used 
some other coordinates, there'd be an extra Jacobian factor. 
We have 

V^T^'^ = e{f)V,TZa + 5{f){d,f)TZa + &U)^,TZ,,ce- (17) 
So, in addition to (^, we have the boundary conditions 



(t^Ai/)rfluid + ^M^surfacc 



= 0. (18) 



Also, when we take the surface tension to be constant: 

□/ (9'^/)(9V)V^5a/ 



^ surface 



97 = -a 6 97, (19) 



where is the trace of the extrinsic curvature of the surface, as seen from outside the fluid (see 
Appendix ^). 

If we have several disconnected surfaces, it is convenient to make the separation / — fi. As 
the surfaces are disconnected, the zero sets of the fi do not intersect. Also, the fi are all positive 
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inside the fluid. Therefore, whenever one of the fi is negative or zero, all the others are positive. 
Luckily, (p^ splits nicely 



i i 

From the form of the gravity solution, we would expect cte/ P to be similar to the thickness of 
the surface. We can estimate it using 

C=- = (20) 

In general, it will be of order and is similar to the surface thickness and ^^fp (if Stt can be 
considered similar to 1). 

For the domain wall of in d = 2 + 1 dimensions, the thickness and surface tension are 
approximately 6 x 2^7^ and cr = 2.0 x ^ respectively. This gives C = which is pretty close to 
the thickness. 

In d = 3 + 1, the domain wall of Q has thickness and surface tension approximately equal 
to 5 X 2^7^ and cr = 1.7 x ^ respectively. This gives ^ = which is also pretty close to the 
thickness. 

For our purposes, it is more convenient to talk about the length scale 

e' = (d+i)e = f = (21) 

Po aic 

2.4 Equilibrium configurations 

In this subsection, we will specialise the general discussion above to the construction of equilibrium 
configurations of fluids with surfaces. We will also derive a simple approach to studying the 
thermodynamic properties of these configurations. 

2.4.1 Solutions for the interior 

We want to find solutions of (||) that are independent of time, which means we need to set (^) to 
zero. This means we need velocity configurations that have zero expansion and shear. In general, 
this would be a combination of a uniform boost and rigid rotation. We can always boost to a 
frame where the centre of rotation is static and the rotation lies in the Cartan directions of the 
rotation group. This gives 

U = -l{dt+^ala), (22) 

where Wq are the angular velocities and la are a set of commuting rotational Killing vectors. The 
important feature is that the velocity is a normalisation factor times a Killing vector (see §2.2 of 

u^'=-iK>', j^K^K^ = -l, V(^X„)=0. (23) 

One can deduce that 

= a^'' = 0, u^'^f,J = 0, = - 

Which leads to 
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One can also show that 





r 






.7. 
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So the velocity configuration ( p^ ) will be an equilibrium solution to the equations of motion 
provided that 

r 

— = T = constant. (24) 

7 

Using the equation of state and (|l|) , this determines all of the intensive thermodynamic quantities 
in the fluid. 

2.4.2 Solutions for surfaces 

The fluid configurations described in the previous subsection have Tji^^jp^^j^^ = 0. Therefore 
This means that (^8|) and ( p^ ) reduce to 

r\f^o = aQ. (25) 



As the pressure is determined by (g4|), this provides a differential equation that determines allowed 
positions of surfaces. Demanding that the surface has no conical singularities turns out to provide 
enough boundary conditions to determine the position of the surface completely (up to discrete 
choices) in terms of the parameters fla, T . 

2.4.3 Thermodynamics of solutions 

We compute the extensive thermodynamic properties of these solutions by integrating the time 
components of the corresponding currents (noting that the current associated with a Killing vector 
Qt- is J'' =T^''Ci^): 



Qx= j dV^x. (26) 

We are assuming that the space-time in consideration is static, so it can be foliated by space-like 
surfaces t — constant with normal dt- In fact, here we will only consider fluids in flat space. 
In particular, also noting that for equilibrium configurations d^f = 0, 



Qc - y Aveu) [{p + rh^K°K-c + pc°] - J dy(5(/) Va7^<. (27) 

Noting that K" = {dt)° = 1 and 1° = 0, this gives 

E - -Qd, - - JdVOif) [{p + Vh^K-dt +V]+ jAVSiD^/dfdjcr, 

La = Qu = J dV0if) [{p + Vh^K-la] , (28) 

S^Qs = jdV6U)[ls]. 

From these quantities, we can compute overall angular velocities and temperature T ther- 
modynamically 

AE = uJadLa + TdS. (29) 

Note that these quantities are different from the local thermodynamic properties of the fluid in 
its rest frame. In particular, the local temperature, T, only knows about the thermal energy of 
the plasma, whereas the overall temperature, T, also knows about its kinetic energy. A priori^ it 
may not seems that these quantities have to be the same as Wa, T from (^ ) and (P^). However, 
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we can show that they are the same by checking that (g9|) holds with uja, T taken from (^2|) and 
(p^). In practice, it is easier to verify the equivalent statement 

d{E^ UJaLa - TS) ^^Ladoja-S dT. (30) 

First, making use of (|^), we see that 

E - ujaLa -TS = -Qk-TQs^- j dVe{f)V + J dVSif) Vdf-dfa. (31) 

Note that the second integral is simply a times the surface area: as we saw in (|l|) the factor of 
y/df-df provides the correct change of measure for the delta function to localise the integral to 
the surface. 

Consider an infinitesimal change of Wa, T. We have 

p + V 

dV — s d{'yT) = dj + 7s dT. 

7 

7-3 d7 = K-dK = K-laduJa- 

From this, we see that (|o|) is satisfied by the contributions from the interior. As the right hand 
side of (|30| ) has no contributions from the surface, we need to check that the surface contributions 
of the variation of (|3^) cancel. 

The change in the surface area can be written as 

dA ~ (f> dAn-w, 



where the integral is performed over the union of the initial and final surfaces, n is a unit normal 
vector pointing into the initial fluid and out of the final fluid and w is some vector field that is 
equal to the outward pointing normal at both the initial and final surfaces. By Gauss' theorem, 
this can be written as 

dA^ I dV\7-w, 



with the integral performed over the region between the two surfaces. The volume element can 
be written as JdV = JdA{n- Ax), with n pointing outwards. As the volume element is already 



infinitesimal, we can replace w with the vector field described in (107), as the difference would be 
infinitesimal, i.e. V-w — Q. Also, as / = on the initial surface, and f + df — on the final 
surface {df refers to the change in / due to the change in ila, T), we have 



df ,^ df 

df 



df.fAx^' + ^d^a + T^dT = 0, 



■Ax = 



VdTdf 
Therefore 

dA = JdVd{f)Qdf. 
So, we can write the surface contribution to the variation of (^ij) as 

d{E ~ LJaLa - T5),urface = -J dVS{f) ^f + J dF<5(/) aO df, 



which vanishes due to (37) 



The thermodynamics of the solution can be summarised by defining a grand partition function 

2gc = Trexp(-^^^). (32) 
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In the thermodynamic hmit, 



We have seen that 



-T In Zgc = E - LOaLa - TS, ^^^^ 

d(TlnZgc) — La du!a + <S'dr. 



ThiZgc= / dVV~ / dAcr (34) 
Jf>Q Jf=o 



and the uja, T are the same as those given by ( p^ ) and (p4|). 
2.4.4 Validity 

We are making several approximations in this paper. First of all, the Navier-Stokes equations are 
merely a long wavelength approximation to the full dynamics of the gauge theory. This will be a 
good approximation provided that the length scale of variation of T and is small compared to 
the thermalisation scale of the fluid (p^. 

Second, we have treated the surface of the plasma as sharp; in reality this surface has a thickness 
of order ^ (see (|20|)). Consequently, our treatment of the surface is valid only when its curvature is 
small compared to 1/^ ~ 1/^' (higher derivative contributions to the surface stress tensor, which 
we have ignored in our treatment, would become important if this were not the case); further we 
must also require that only a small fraction of the fluid should reside in surfaces. This boils down 
to demanding that all sizes are much larger than £f . 

Thirdly, we have ignored the fact that the surface tension is a function of the fluid temperature 
at the surface, and simply set a = cr{Tc)- This is valid provided that T /Tc « 1 at all surfaces. 
When this is the case, the pressure will be small compared to pc (see (||)). Then, ( ^7|) tells us that 
the extrinsic curvature of the surface must be small compared to 1/^, which is the same as the 
previous condition. 

We can estimate the scale over which thermodynamic quantities vary as the distance over 
which the fractional change in the temperature is one. As the temperature is proportional to 7, 
we should demand (schematically) 

1 l-v^ , 

-j] 77- ^ tmfp • 



||Vln7|| 

At temperatures close to Tc^ where our other approximations are valid, we have Zmfp ^ There- 
fore, we require that that 

This will be true if the speed of the fluid is much less than the speed of light and the angular 
velocities are much less that 1/^'. 

2.4.5 Dimensionless variables 

It is convenient to rescale the variables as follows. First we rescale all lengths and times by ^' (i.e. 
working in units where ^' = 1) 

X = ^'i, = ^ . (35) 

We measure temperature in units of Tc and further rescale extensive thermodynamic quantities 
by aTc". 

Then, (|l]) becomes 

_ 1 e, (37) 



sur 



face 
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and ( ^ ) becomes 



Jf>o L J ^/=o 

-T In Zgc = - 5,La - fS, (^^^ 
d(TlnZgc) = ia duja + S dT. 

From now on, we will suppress all tildes and work entirely with the new variables. 
2.4.6 Summary 

We can summarise the construction of equilibrium solutions as follows: 
The fluid velocity is given by 

M*' = 7if ^, where K = dt + uJads , 

(39) 

with (j)a being a set of angular coordinates such that (j)a — > 0a+Ca are a set of commuting isometrics 
and uja the angular velocities. 

The thermodynamic properties of the fluid are specified by 

T = iT. (40) 

All other local thermodynamic properties can be computed from the equation of state (|^) and the 
relations (|^). 

The position of the surface is specified by a function / that is positive inside the fluid and 
negative outside. It is determined by 

[(7r)'^+^-i]_f_-e, (41) 



with 8 given by (109), and the condition that the surface is closed without conical singularities. 
The overall thermodynamic properties of the solution can be computed from 



TlnZge^/ dV[{^TY+^-l]- j dA, 

^/>o ^/^o ^^^^ 



//>o 

-T In Zgc = E - LOaLa - TS, 

d(T In Zgc) — La duja + S dT. 



3 Black objects in arbitrary dimensions 

As mentioned in the introduction we study the horizon topologies of black objects in SSAdSd+2 
through the dual fluid configurations which solve the d dimensional relativistic Navier-Stokes 
equation. In |^ exact disc (-B^) like plasma configurations were obtained in 2+1 dimensions. In 
one higher dimension (i.e. in 3+1 dimensions) we can expect the existence of the solution x R^, 
which is topologically a cylinder. Now we can bend the cylinder into a ring. The fact that such 
a ring solution exists were shown (numerically) in Q and its thermodynamics were explored in 
. This ring solution can be constructed perturbatively from the cylinder (a suitable expansion 
parameter being the ratio e = ^, R being the radius of the cylinder and C being the distance of 
the cylinder from the origin). The leading order solution, called the thin ring, (correct up to 0{e)) 
was obtained analytically in |^ . 

This method can be used to construct rings and other exotic plasma configurations in higher 
dimensions. For example in one more dimension (i.e. 4+1 dimensions) we can have a cylinder 
solution with topology x {B^ being the ball type solution in 3 + 1 dimension) then we can 
bend this cylinder to form a ring. In this way we would be able to obtain ring type solutions 
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from the plasmaball configuration in one lower dimension. Further besides the ring we can also 
construct other configurations by a similar method. Again for concreteness let us consider the 
example of 4 + 1 dimension where besides the x topology we can also have a topology 
X constructed out of the B^ solution in 2 + 1 dimensions (two dimensions lower). Now 
we can bend the B^ x configuration along two directions in the R^ to form B^ x x S^. 
Another way to think of it is that in 3 + 1 dimensions the B^ x topology existed. So (just like 
the construction of the cylinder) in 4 + 1 dimension we have a topology x 5^ x R^. Now we 
can bend this cylinder along the R^ to form the topology B^ x x . Thus in d space-time 
dimensions we can use this method to study a topology 

^{d-l-n) ^ j,n ^ Q(d-l-n) X S"" . (43) 

n times 

Note that there can be solutions with other topologies in d dimensions which we will not be able 
to capture by our method, therefore we shall look at these topologies only. Here n cannot be 
arbitrary. Naively we would expect n < d — 3 so that the ball is at least a B^ |^. However there is a 
stronger bound on n. We recall the fact that for the ring (or more exotic objects as above) to exist 
we must have non-zero angular momentum in the plane of the 5*^ (s) . However in d dimensions we 
can independently turn on angular momentum only along (^^) directions for odd (even) d. 
This is because the group of spatial rotations in d dimensions is SO{d — 1) which has rank 
(^^) for odd (even) d. Thus n should not be more than the number of Cartans of the (spatial) 
rotation group in a particular dimension. Hence n should satisfy 

n = 0, for d = 3. 

n < — ^ — , for odd d greater than 3. ^^^-j 

d-2 ^ 

n < , for even d. 

- 2 ' 

Here note that d = 3 is a special case because for this dimension rt < d — 3 is a stronger bound 
than n < Also note that a new topology of the plasma configuration is obtained for every 

odd d (> 5). Thus in an even d we only have the solutions that existed in d — 1 dimensions. 

Here we shall analytically show (in an approximation analogous to the thin ring approximation) 
that these topologies exist as plasma configurations which are solutions of the relativistic Navier- 
Stokes equations expressed in the form (^). The most suitable coordinate system for studying the 
general topology ( ^ ) is the one in which we choose {tq, (j)a} with a = 1, ... n as the coordinates on 
the n planes containing the n S^s and in the rest of the space (which exists exists when d > 2n-|- 1) 
we choose spherical polar coordinates {r, 6'i, . . . , 0{d-2n-2)}- The coordinates (j)a and 9j are angular 
coordinates; the coordinates {6j} with j — 1, . . . , (d — 2n — 2) may be taken to be the coordinates 
on a unit sphere in a d — 2n — 1 dimensional space. The angles 9j will be absent if d = 2n -I- 2. 
The range of all the radial coordinates {r,ra} are as usual [0, oo) while the angles (j)a G [0,2tt). 
Besides these spatial coordinate we denote the time coordinate by t. 

The metric is given by 

n n 

ds' = -dt' + J2 ^'^'t'l + E d*^' + + '^^'"'""'^df2^,_2n-2) {0l,02, • • • , ed-2n-2) (45) 

a— 1 a— 1 

where dr2^^_2n-2) round metric on a unit sphere in d — 2n — 1 dimensions. In such a space 

we consider the fluid surface to be given by 

/ = /i(ri,r2,...,r„)-r = 0. (46) 



If the ball is a then we would be considering hollow solutions; however we are unable to capture the 
description of those solutions by our method. In addition, such hollow solutions might not exist hi all dimension. 
In 2 -|- 1 dimensions the existence of such hollow solutions (the annular ring) was reported in [pj; while there is 
strong evidence (see Isl B|) to suggest that such hollow solutions do not exist in 3 + 1 dimensions. 
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For n = we have a special case because then there is no r coordinate. In this case we 
consider / = h(r2, • ■ ■ , rn) — ri = as the equation of the surface. 

The velocity is given by = 7(1, wi, . . . , cj„, 0, 0, . . . , 0), with 7 = (1 - Y^'^^^ujlrD^ being 
the normalization factor. Note that here we are not considering any angular velocity in the 9i 
directions. The topologies that we explore could also be spinning along the 6i directions. However, 
if that were the case the zeroth order solution would not be simply a round ball times a plane and 
we would lose analytic control. It might be possible to turn on infinitesimal angular velocities in 
these directions, but we will not consider that here. 

The equation for h{ra) that follows from ( |37| ) is given by 



/ n \ n 



(47) 



„b=l 



where m — d ~ 2n — 2. In order to obtain (^ ) from ( |37|) we set the typical length scale of 
the problem ^' to 1 by suitable choice of units and we are measuring temperature in units of Tc, 
as described in §2.4.5. For n = the sum starts running from 2 (instead of 1) in the above 
equation. Now we separately consider the space of {r,ra} (only the first 'quadrant' of this space 

is physical because {r,ra} G [0,oo)). We now shift the origin to a new point ~E such that it is 
given by 

1 = {0,CPl,CP2,...,CPn) (48) 

PaS are projectors along various directions, so that ^ Pa ^ ^ a-nd C is the magnitude of the 
vector ~C . Again in the special case of n = we take ~^ = {CPi, CP2, . . . , CPn)- Let {xa} be 
the new shifted coordinates such that 

ra=CPa+Xa. (49) 

As is apparent the coordinate r (if it exists) remains unchanged by this coordinate change. 
Now we perform the following scaling 

UJa^eWa, C = —; h = y{{xa}). (50) 
e 

Then ( ^t]) at leading order in e (which is e'^) reduces to the equation 



a=l 



(51) 



\ a=l / a,b=l 



Now (^ is satisfied by the function 

yi{xa})=(^R'-j2^l^ ' (52) 
provided the following equation is satisfied by the parameters 

T^^^-( ^''%^^ + ^ ) (l-E(^oW) ^ (53) 
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Also the equation ( ^ ) at 0(e) yields 



(d+l)r'^+lEa(<^a4Xa) 1 ST ^ ^ Q (54) 

In the above equation if we substitute (^2|) and set the coefficients of Xa to zero then we get (after 
using (H)) 

" (^oPa)'(((d-n~2) + i?)(d+l)+n)' 

Note that here we have n + 1 equations for n + 2 parameters R, io and n PqS. However there is 
one constraint among the PaS (namely J^Pa ~ 1) which gives us the correct number of equations 
for the parameters to be determined. 

Again for the special case of n = the chief results (^2|), ( ^3| ) and ( |55| ) remain unchanged. 
However we have h = CPi + y{{xi}). As a result the 0th and 1st order equations, ( pi] ) and (|5^), 
are changed respectively to 



(56) 



a,h 



and 



(<i + l)T^+iE„WP„f„i„) 



As mentioned before here the sums run from 2 to n. Also in this case besides equating the 
coefficients of XaS to zero in ( |57| ) we also have set the coefficient of {B? — Ea-^a)^ to zero to 
obtain (p5| ). 

In (pq) we find a very curious fact about the speed of the class of solutions that we analyze 
here. The velocities WaPa^o reach a maximum when — >■ and the maximum value is given by 

wT'^PJa = ^ (58) 

This value is consistent with the maximum speed for the ring in 3 + 1 dimensions quoted in 
Also note that for the ring [n — \) and large d this maximum value goes as ^; this is unlike 
(although consistent with) the behavior of the ring in asymptotically flat space where this goes 
as (see Q). However, this limiting value occurs when i? — 0, where our approximation 
of the surface as having no thickness breaks down. Nevertheless, for large space time dimension 
the behavior of black holes in asymptotically flat spaces and that in asymptotically AdS spaces 
are expected to be similar (see [^). In the light of this fact we may conclude that our fluid 
approximation is unable to capture this phenomenon correctly, unless there actually exist a better 
bound in flat space (which would go as ^ for large d). 

In the coordinates that we have used in this section, the first derivative of our solution 



dby{xa) ~ — Yl ^ ^ ) near J2^a — P^ ■ In fact (from our analysis so far) it is 

unclear whether a consistent perturbation theory can be performed in the e parameter about the 
solution (^2|). In the later sections we shall move to better coordinates and exhibit the existence 
of a well controlled perturbation theory about (^2|). In certain special cases we shall explicitly 
compute the ffist correction to ( ^ which occurs at O(e^). 
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Summary^ 

Here we have constructed a class of fluid conflgurations to the d dimensional Navier-Stokes equation 
in the generalized thin ring limit. To leading order in the parameter e these fluid conflgurations 
are given by, 

^id-l-n) X j,n ^ ^(d-l-n) x 5^ X , 

v ' 

n times 

where n < [^^]- These conflgurations are parameterized by the radius of the ball (i?) and the radii 
of the various S'^s (ioPa). In the generalized thin ring limit locally these configurations are like 
filled cylinders with the topology B^"*-!-") x M". Then we can bend the different directions in M" 
into S^s in a controlled way with a perturbation expansion in e. Now the intrinsic fluid parameters 
(namely the temperature (T) and the angular velocities (wa)) are related to the parameters of the 
fluid conflguration {R and ioPa) by the force balance conditions. The pressure along the radial 
direction of the ball is balanced by the surface tension. This condition yields 

On the other hand the pressure along the radial direction of the S'^s is balanced by the centrifugal 
force. In order to obtain this force balance we require these conflgurations to be rotating (at least) 
in the planes in which the S^s lie. For the sake of simplicity we have turned off angular velocity 
along any other direction. This force balance determines the angular velocities to be 



" {£oPafm-n-2) + R){d+l)+ny 

These fluid conflgurations are dual to the horizon topology S'('*^") x T". Thus by exploiting the 
AdS/CFT correspondence we indirectly conflrm the existence of such exotic horizon topologies. It 
is possible to generate a perturbation expansion (in the e parameter) about these solutions as we 
shall demonstrate in some explicit examples in the later sections. Also the thermodynamic prop- 
erties of the fluid configurations directly map to that of these exotic black objects. This provides 
us with a opportunity to study the thermodynamics of these black objects without performing a 
direct gravity calculation. 



4 Rings 

In this section we shall analyze the topology B'^~^ x S^. This topology is the special case of ( ^ ) 
(for n — 1). We shall study this ring type solutions in 3 + 1 and 4+1 dimensions in detail. 
We use regular coordinates to set up a well controlled perturbation expansion in the parameter 
e to compute corrections to the thin ring. The thin ring solution in 3 + 1 dimensions has been 
well studied (including the thermodynamic properties) in |^ . Here we present the first correction 
to that solution demonstrating the fact that it is possible to construct such rings as a series in 
the e parameter. This method can in principle be generalised to higher dimensions (we present 



the results in 4 + 1 dimensions in p.2| ) . All these solutions are consistent with the general case 
discussed in sj^ to zeroth order in e. 

4.1 Rings in 3+1 dimensions 



As pointed out previously in 
are those that are used in [§l 



3^ the coordinates that we used for our general discussion (which 
I) are not suitable for the perturbation theory. Therefore, for the 



construction of a regular well controlled perturbation theory, we move to some regular coordinates. 



^While this work waa.^being completed we came to know that similar results were also obtained by Camps, 
Emparan and Haddad 
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Here we use the coordinates {t, p, 0,(j)}, t being the time coordinate, which are related to those 
used previously by 

r — £ + p cos 9, z = psm9. (59) 

The metric is given by 



'dr + dp-" + p^dO^ + (£ + pcos 



\2 J J,2 



(60) 



Here £ is a number which we shall determine in terms of the parameters of our system order 
by order in perturbation theory. Physically C is the distance of the center of the cylinder (our 
0th order solution) from the origin. We simply redefine this center to be our origin in the above 
metric. 

These coordinates are described in figj2. 




Figure 2: Cross section of the 3+1 dimensional ring. The curved arrow labelled (j) indicates a 
direction that has been suppressed. 



The velocity vector is given by — 7(1, 0, 0,a;), where again 7 = (l — w^(£ + pcos(6'))^) 
is the normalization constant. In this case we consider the fluid surface to be given by 

f = g{0)-p = O. 

Then ( ^7|) reduces to the following differential equation for the function g{6) 



1 



(1 - Lo^C + cosi9)g{9)rf' {C + cos(%(0)) (5(^)2 + g'i0)^f' 

+{2C + 3 cosi9)g{9))g'{9f + g{9f s\n{9)g'{9) + g(9) {g{9){L + 2 cos{9)g{9)) 
-iC + cos{9)g{9))g"i9)))^0. 

Now we perform the following scaling: 



sm{9)g'i9f 



(61) 



oj — e w 



g{9) = R + eg,{9) + e^g2{9) + e^g3{9) 

£ = -io + h + + 
e 



(62) 



where e is the small parameter with which we wish to perform the perturbation. Here gi , 32, 53 are 
functions to be determined and ^0j^i;^2,^3 are to be expressed in terms of the fluid parameters. 
Also note that to zeroth order g{9) = R, which is the same solution that has been obtained in 
Also the solution in (as in the general discussion above) was true up to first order in e. This 
implies that the first order correction to g{9) (i.e. gi{9)) should vanish, as we shall shortly show. 
To first order in e (|6l| ) reduces to 



1 

R 



rp5 



0. 



(63) 
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To higher order in e we obtain differential equations for (?i,g2,ff3 etc. These differential equa- 
tions are of the general form 

g.m+g':{e)^S,{e). (64) 

where S{9) is the source which is determined at a particular order once the the complete solution 
up to one lower order is completely known. Also note that the homogeneous part of the equation 
is the same at all orders. 

The equation that we obtain at first order is 

Solving the above equation we obtain 

31(6*) =Ci cos(6l) +C2 sin 61 

i?2 (((1 - £2^2^ _ 5r5£2^^2\ (^^g^^) _^ 26'sin(6')) - 20THliiwA (66) 



where the Ci and Ci are the integration constants to be determined by the boundary conditions. 
We must also remember that the constants R, io, ii are also to be determined in terms of the other 
fluid parameters, viz. K and w. 

Now from physical considerations we shall demand that the surface should be a closed surface. 
This results in the boundary condition 

g[{0) = Q; 5'i(7r) = 0. (67) 

The first condition demands C2 = 0, while the second condition yields 

5T^£gi?w2 



(1 - £2u;2)7/2 



1 = 0. (68) 



Here the relations (|63| ) and ( |68| ) may be used to express R and £q in terms of the fluid parameters 
T and w. However for performing the calculations it is more convenient to do the reverse. We 
then have 



Rj \6 + 5R J ' £oV6 + 5i? 

Plugging back these relations into ( |66| ) we find g{9) up to order e, which is given by 

g{9)^R + e(cicos{9)-^) +0{e^), (70) 



where Ci and £1 are still to be determined. 

Now once we start including corrections we should consider a redundancy in description of the 
ring which we have to remove by proper gauge fixing. This pertains to the fact that we haven't 
defined C properly yet. Vaguely, it is r coordinate of the center of the ring. This becomes ill- 
defined once we take into consideration the corrections to the thin ring, which was a circle in the 
p-9 plane to zeroth order. This is taken care by the following gauge-fixing condition 

gi9)cos{9) = 0. (71) 

In words, this condition states that the average r coordinate of the surface in the r-z plane is C. 
This condition implies Ci ~ 0. The constant £1 will only be determined at next order in epsilon 
just as £0 was determined at 0{e). Therefore we now proceed to the calculation of the second 
order corrections to e. 
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The equation for 52 (^) is given by the coefficient of in (|6l|). After plugging in the solution 
for gi{6) we obtain 



+Rcos{e){2£i{5R + 12) + R{lbR + 22) cos(0))) . 
The solution to the above equation is given by 



(72) 



52(6') =C3 cos(e') + C4 sin 6* 
R 



■ 60e{R + 1) + 19) cos(^?) + R{15R + 22) cos(20) (73) 

-eiiibR + 12)0sin(0)) - 3 ((4 - 10R)£l + 2044(i? + 1) + R'^{WR + 22))) 

where again C3 and C4 are integration constants to be determined. The boundary conditions (i.e. 
52(0) = = 52(''')) again imply C4 = and the condition 

hnR^5R + 12) _ 

mliR + 1) ^ ' 

The above condition implies = 0. At this point the e order solution is completely determined; 
and we find that corrections to g{9) and C all vanish. However we go ahead further to compute 
the O(e^) correction completely as it will provide us with the leading order correction. Using the 
fact £1=0 and C4 = we find g{9) to be 



g{9) = R+C3Cosie)e^ 

R (i?2(i5i? + 22) cos(26i) - 3 ((15i? + 22)R^ + 204^2(i? +1)))) 2 _^ ^3 



eoeliR + 1) 

(75) 

Now the condition (Q) again implies C3 = 0. Again in order to determine £2 we have to perform 
one higher order calculation. We can then determine £2 by imposing the boundary conditions on 
93(9) ■ Here we intend to present only the leading order corrections and therefore do not specify 
the details of the third order calculation. However the value of £2 that we obtain is given by 

_ i?2 (225i?2 + 380i? +92) 
" 40£o(5i?2 + 17i?+12) ■ 

Note that the denominator in the above expression never vanishes for positive values of R. 
In summary we can write 



n i?^(2(5i?+12)(15i? + 22)cos(26')-15(3i?(25i?+64) + 124))e2 ,3. 

3^^^-^+ 120£^,{R+1){5R+12) + ^ > 

1 R"^ (225R^ + S80R + 92) 
^ = 7^° + - 40M5i.^ + 17i^+12) 



(77) 



As mentioned earlier R and £q are related to the fluid parameters T and w through the inverse 
of the relations (|69|). We can carry forward to arbitrary order in e in a perfectly well controlled 
fashion. 

Finally we would like to mention that it is possible to obtain corrections to the radius R from 
([ttI). However, just like the center, the notion of the radius has to be redefined for the corrected 
solution. This we may do by defining an average corrected radius, 

i?avg = / g{9)d9. (78) 
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Using (|77|) we can compute i?avg to be 



Ravg — i? — 



(i?3(3i?(25i? + 64) + 124))e2 
8(^§(i? + l)(5i?+12)) 



We present a plot of this corrected solution in fig.|2 





9(6) 




circle 




(radiuEi = 











(79) 



Figure 3: A plot of the corrected solution ([T^) for the 3+1 dimensional ring with — = 14 and 
i? = 7. 



The thermodynamic properties of the solution can be computed from ( pq ) as follows 

raie) 



T In Zgc = 2tt 



We find 



/•T ray.") 
de dpp{C + pcos0)iK-f^ -1) 
Jo Jo 



deVgW+TWi^ + 9(0) cos ( 



o — 

e 

27r2i?i6/V5i? + 6(i?(5i?(105i? + 353) + 2018) + 792) 
8%/54(i? + l)"/i"(5i?+12) 



O 



(80) 



_ 27r%i? 27r^fl3(15i? + 22) 
Tln^g,- _ + _____, + 0(, ), (81) 

The other thermodynamic properties can be found by differentiating this (|33|): 

27:^£oRi5R + 7) 27:^ R^{R{5R{105R + 353) + 2018) + 792) . 
" ^ 8(4(i? + l)(5i?+12)) ^ ' 



(82) 



_ 2tt^£IRV5R + 6 27r2i?V5i? + 6(15i? + 22) . 
+ 40(i?+l) 

Note that to leading order these expressions for energy, entropy and angular momentum match 
with that presented in |^ after performing the necessary variable transformations. 
We present a plot of the thermodynamic properties of this ring in fig.^ 

4.2 Rings in 4+1 dimensions 

We now analyze the ring in one higher dimension i.e. 4 + 1 dimension. The construction is exactly 
parallel to that in 3 + 1 dimensions and therefore we skip most of the details and specify only the 
result. 
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Figure 4: Plot of entropy, S, vs. angular momentum, L, at fixed energy, E = 10,000, for the ring 
in 3+1 dimensions. The dotted portion of the curve represents the region where R < 1 and hence 
lie outside the surface tension approximation. 



The coordinates that we use here are {t, p, 9, (f>i, 4)2}, which are related to the old coordinates 

ri=C + pcos9, r2=psm0, (83) 



by 

with the metric 



ds^ = -dt^ + dp^ + p^de^ + {C + p cos 9)^d(j)j + [p sin ( 
These coordinates are described in figj|. 




ri 

<t>i 



(84) 



Figure 5: Cross section of the 4+1 dimensional ring. The curved arrows labelled (j) indicate a 
direction that has been suppressed. 

The velocity is given by u'^ — 7(1, 0, 0, 0) where the normalization is given by 7 = (1 — 
uj^{C + pcos(6'))2)-5. We take the fluid surface to be / = g{6) - p = 0. Then the equation (|^) 
reduces to 



2^6 



1 



3/2 



'(9)2 



3/2 



(gior + g'ief) {g{e){2£ + 3 cos(%(0))- (^s) 



csc(6l)(£cos(6l) + cos{29)g{9))g'{0) - (£ + cos{0)g{0))g'\0)) = 0, 
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Now we plug in the scaling with e as in the previously discussed cases: 



9{e)=R + egM+(^92{0) + 



(86) 



e 



e^2 + . . . 



Then, performing an analysis exactly the same as the one performed in §[4.1|, we obtain 



m = R+ 



(5 (81i?2 + 432i? + 572) cos(26i) - 3 (243i?2 + I040i? + 1108)) ^ , 3. 

240^2 (9^2 + 44^ + 52) ^ +^^(0 



1 i?2/'8j^^2_,_24o_R+ii6) 
^ = 7^"+ 30M9i^^+44i?+52)^ + ^(-) 



(87) 



with R and being expressed in terms of the fluid parameters T and w (implicitly) by the 
following relations 

216(i? + 2)4 _ 1 



i?(6i?+ 13)3 ' 



4V6i?+ 13 

Even in this case we note that all the 0{e) corrections vanish. 

Finally the average radius of the curve in (|87| ) (-Ravg as defined in ([zsj)) is given in this case by 

(i?3(243i?+554))e2 



80 (^^(9i? + 26)) 
We present a plot of this corrected solution in fig.g. 



(89) 
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Figure 6: A plot of the corrected solution ( |8^ ) for the 4+1 dimensional ring with ^ 14 and 
i? = 7. 



The thermodynamic properties of the solution can be computed from ( pq ) as follows 



rinZgc = (27r)^ 



d0 dpp'^sm9{C + pcos9){K-f'^ -1) 
Jo Jo 



Ae^g{ef + g'{eY{c + g{e) cos e)g{e) sine 



(90) 
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We find 



TlnZf^c = — ( -6i? + \/lTT2V36R + 78 - 18 



+ 



-3216V36i? + 78i?^/^ - 243\/36i? + 78i?^^/^ 



(91) 



1354(-R + 2)5/3(9i? + 26) 
-1524V36i?+ 78i?^^/^ + 1458(i? + 2f/^R^ + 9630(i? + 2)^/^7?^ 

+21888(i? + 2f/^R - 2288V36i?+78v^ + 17160(i? + 2)^/3^ e 

For expressions of the energy, angular momentum and entropy refer to Appendix We present 
a plot of the thermodynamic properties of this ring in fig.^ 
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Figure 7: Plot of entropy, S, vs. angular momentum, L, at fixed energy, E — 1,000,000, for the 
ring in 4+1 dimensions. 



5 'Torus' in 4+1 dimension 

Here we analyze the solution with the topology x x S-^ which for the lack of terminology we 
refer to as the 'torus'. Although the perturbation theory for the torus is almost exactly parallel to 
that for the ring, however there are certain differences as far as imposing the boundary condition 
is concerned. We consider spatial part of the 4+1 dimensional space to consist of 2 independent 
planes. We turn on two independent angular velocities (a;i,W2) along a direction orthogonal to 
these planes. The two S^s of the topology x x lie on these two planes. Initially we put 
polar coordinates and r2,4>2 on these two planes. Further in the (ri,r2) plane we shift to 

coordinates p, after a shift in the origin by the vector = {C cos(x), sin(x)) (expressed in the 
(ri, coordinates). The various coordinates have been represented in figjsj Thus finally we work 
with the coordinates {t, p, 6*, 02}, which are related to the old coordinates by 

ri = £cosx + pcos6', r2 = >Csinx + psin6', (92) 

with the metric 

ds^ = -dt^ + + p^de"^ + (£ cos X + /o cos OfA(l)l + (£ sin x + p sin efd(j>l (93) 

Here the velocity is given by — 7(1, 0, 0, wi, 012) where again the normalization is given by 
7 = (1 — (£cosx + pcos9)'^ujf — (£sinx + /9sin6')^aj|)^2 . We consider the fiuid surface to be 
given by / = g{0) - p = 0. 
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Figure 8: Cross section of the 4+1 dimensional torus, 
direction that has been suppressed. 



The curved arrows labelled 6 indicate a 



Now the differential equation satisfied by g{9) (which is (p7|) for the present case) is given by 



2^6 



(1 - (£ cos X + 9[e) cos ey ~ uj^{C sin x + 9{d) sin 9)^) 



- 1 - 



{9iO)+9"{d))9'{0? 



(<7(0)2+5'(0)2) f cos{9)9'{9)-9{0)sm{9) cos(%(^^) + sin(%'(0) 



3/2 



g(0)2(£W +1^'/' \sm{9)9{9y + Csin{x)9iO) cos(%(0)2 + £ cos(x)5 W 

9"i9)-9iey 



(94) 



9ioy 



= 0, 



(95) 



Now we again consider the following scaling: 

LUi = e wi 

9{9)^R + egy9)+y 92{9) + ... 
e 

X = Xo + e Xi + X2 + ■ • ■ 

We shall determine the unknown functions in a similar way as we did for the ring. However there 
is a crucial difference between the two. Firstly here we have one more parameter (since we have 
two angular velocities instead of one) . Secondly the boundary condition that we have to impose on 
9{9) is different from the previous case because here we are dealing with a different closed surface. 
Although physically it is the same criterion - the fact that we should have a closed surface, the 
mathematical formulation of the statement is different as we shall now describe. 
Instead of the boundary condition (|67|) we should use the condition 



5(0) =5(2^); ff'(0)=5'(2^)- 



(96) 



which is the statement that we should have a closed curve in the p-9 plane and that the curve must 
close in a regular fashion such that the derivatives on either side of the point of closing (which we 
take to be = 0) are equal. As the differential equation is second order and periodic, this ensures 
that all higher derivatives are continuous at 9 = and the solution is fully periodic. 
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Besides the boundary conditions we will also have to fix the ambiguity regarding the center of 
the torus just as we did for the ring. However unlike the ring the center here does not lie on the 
r2 axis. Therefore in order to fix the center we will have to use two conditions namely 

27r /•27r 

g(e) cos{e) =0, / g{0) sm{9) = 0. (97) 



In words, these conditions states that the average ri coordinate of the surface in the ri-r2 plane 
is >Ccosx and the average r2 coordinate is Csinx- Then proceeding in the same way as the ring 
(after including the above modifications) we find the following result: 

9(0) = R- 72 {e{R+l)) + + + ^^^^ + " 

+4(9i? + 17) cos(2(^? + xo))) csc2(2xo)) + O (e^) 

4 i?2(81i?2 + 150i? + 8cos(4xo) + 57)csc2(2xo) ^,3, (98) 

"^ = 7 + 48Mi^ + i) ^ + 

i?' cot(2xo) 2 

Note that just as for the ring all the order e corrections vanish. Here R, £q and xo are again given 
implicitly in terms of the fluid parameters T, wi and 'W2 by the following relations 



27(i?+l)^ 1 / 1 1 / 1 .9g^ 

i?(3i? + 4)3' ^0 cos(xo) V 8 + 6i? ' 4sin(xo) V 8 + 6i? ■ 



Note that this entirely matches the general results (^3|) and (|55|) for d = 5 and n = 2. 

Again in this case, the average radius of the closed curve ( P8|) (with i?avg as defined in (|7^)) is 
given by 

R -R ^MfiM27i?^ + 74fl + 5l)csc^(2xo)) , ,..3. 

^avg - i? 8(fg(i?+l)) + ^ 

We present a plot of this corrected solution in fig.^. 

We expect our construction of the torus solution to break down when £0 ^ R. Which is 
reflected in the fact that ^"^^ ^'^^ ^ proportional to . 

The thermodynamic properties of the solution can be computed from (pq) as follows 



rinZgc = (27r)^ 



d9 dpp{Ccosx + pcos9){Csmx + psm9){Kj'^ ~1) 
Jo Jo 

- / dOy/ g(9)'^ + g'{9Y(C cos x + g{d) cos 9){C sinx + g{9) sin^ 



(101) 



We find 

27r3f2^sin(2xo) 



ThiZ„c = 



(-3i? + l%/9i? + 12v^ - 6 



3e2 

+ w"^ (^(-^ + l)'^'(^(3i?(27i? + 98) + 401) + 208) (1^2) 

- ^(3i? + 4)V9i?+ 12(i?(27i? + 62) + 39)) + O (e^) , 

Again the expressions for the energy, angular momentum and entropy are given in Appendix 
I We present a plot of the thermodynamic properties of this torus in fig.p^ 
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20 40 60 80 100 120 140 



Figure 9: A plot of the corrected solution ( p8|) for the 4+1 dimensional torus with — = 100, 
Xo = f and i? = 10. 




Figure 10: Plot of entropy, S, vs. angular momentum, Li, at fixed energy, E = 1,000,000, and 
different values of for the torus in 4+1 dimensions. Here again the dotted portion of the curve 
represents the region R < 1 which is outside the validity of the surface tension approximation. 
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6 Discussion 



In this paper we have explored an interesting prediction of the AdS/CFT correspondence which 
connects the dynamics of classical gravity to that of fluid dynamics in appropriate limits. This 
duality has been exploited by several authors to make interesting gravitational predictions about 
the properties of the dual fluid (see |12 for a review). In it was shown that for certain 
special cases it is possible to reduce Einstein's equations on the bulk side to the Navier-Stokes 
equations governing the boundary fluid dynamics. In this process they have established a one-to- 
one correspondence solutions of fluid mechanics and gravitational solutions. 

Here we have used this correspondence in a somewhat reverse fashion to prove the existence of 
a class of very exotic horizon topologies in higher dimensional AdS-like spaces using solutions of 
the equations of fluid mechanics. However, perhaps due the lack of the corresponding gravitational 
solutions, we do not have such a precise map between the horizon and the boundary. Nevertheless 
it is an interesting question worth pursuing in future. 

Also we would like to emphasize one interesting observation regarding the velocities of rotation 
of the solutions that we have studied. In (^8|) we pointed out that maximum velocity of rotation 
along any of the planes has a upper bound independent of the energy (although the i? — >■ limit in 
which this bound was obtained lies outside the validity of the surface tension approximation) . For 
the case of the ring (for which n— 1) this upper bound is given by {{d — 3){d + 1) + 1) ^ (which 
for d = 4 yields consistent with that obtained in 0). Note that for large d this upper bound 

behaves as ^. We may compare this with the corresponding result in flat space presented in 
where the bound is (d — 1)^3 ^ -1= for large d. 

In the authors had proposed that the black objects in SSAdS spaces with large number 
of spatial dimensions may qualitatively behave the same way as the corresponding black objects 
in flat space. This is could possibly be understood as follows: The SSAdS space has one special 
direction, namely the radial direction, in which Hawking radiation from the plasmaball bounces 
back. In the other directions, the Hawking radiation escapes to infinity, much like in flat space. 
As a result, plasmaballs show behaviour that is somehow in between flat space black holes and 
large AdS black holes. For example, they have a negative heat capacity and eventually evaporate 
completely. As we increase the number of dimensions, the flat directions outnumber the AdS-like 
radial direction more and more. Hence, one might expect the black object living in SSAdS spaces 
to behave more like their flat space counterparts as one increases the number of dimensions. For 
example, at fixed energy, plasmarings have an upper bound on the angular momentum in five bulk 
dimensions, like in AdS jl^, but we see no such bound in higher dimensions, like in flat space. 

With this in view, it might be that there existed a bound for the black objects in SSAdS spaces 
that went as for large d and our fluid analysis gives us an incorrect value of this bound (as our 
approximations breakdown in the limit in which the bound is obtained) . Otherwise it might also 
be that the qualitative similarity of the fact that a universal upper bound actually exists in both 
cases is all that we should expect, there being a difference in the exact value of this bound. 

It might also be interesting to study the time dependent fluctuations on these solutions to 
understand their stability properties. Like many fluid configurations these might suffer from 
Plateau-Rayleigh instability that would indicate Gregory-Lafiamnie type instability on the grav- 
ity side |l^ (see po|-p2[ for recent work connecting the Plateau-Rayleigh and Gregory- 
Laflamme instabilities). Since Plateau-Rayleigh type instability is very commonly found in narrow 
strips of moving fluid (like the one that we consider here) it may be true that with our methods 
we may never be able to capture stable conflgurations. 

In five dimensional gravity there are a number of interesting solutions with disconnected hori- 
zons, such as black Saturns with a black hole surrounded by a black ring, and di-rings [p4|-p7|. 
These solutions play a significant role in the phase diagram [^. We note that the construction 
of such solutions in fluid mechanics is trivial - disconnected lumps of plasma do not affect each 
other at all. This may sound surprising from the gravitational perspective, as one would expect 
black holes to attract each other. However, in the AdS soliton background the spectrum of the 
graviton has a gap Agap ^ S,^^- Therefore one would expect the gravitational attraction to die off 
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exponentially rather than the usual power law, e~^sap»' ^ e~^^^. As the surface tension approxi- 
mations is only valid when the separation of these lumps of plasma is much larger than the surface 
thickness, r ^ ^, the gravitational attraction would be negligible. 

One issue concerning Saturns is thermal equilibrium between the components . This would 
require the temperature and angular velocity of the components to be equal (see also ) . It was 
pointed out in that such solutions do not exist in 2+1 dimensions in this framework. However, 
demanding thermal equilibrium between the disconnected components is only reasonable when the 
time scale of energy /angular momentum exchange between the components is much shorter than 
the rate of Hawking leakage to infinity. This limit may not be compatible with the demand that 
the separation of the components is much larger than the surface thickness. 
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Appendices 

A Extrinsic curvature 

Suppose we have a timelike surface with unit normal vector n pointing toward us (spacelike surfaces 
will require some sign differences). The induced metric on the surface is 



(103) 



The extrinsic curvature is given by 



(104) 



We have to be a little careful with the last expression. It agrees with the first expression when 
projected tangent to the surface. The first expression has vanishing components normal to the 
surface. The normal components of the second expression depend on how we extend n off the 
surface. 

The conventional choice for extending n is as follows: at each point on the surface, construct 
the geodesic that passes through that point tangent to n and parallel transport n along it. In 
other words 

nf^V^n" = 0. (105) 



This ensures that the second expression in ( |104| ) has vanishing components normal to the surface. 
The other normal component, n^V^n^, vanishes due to the normalisation of n. 

For the surfaces given by f{x) = 0, considered in §2.3, the unit normal on the surface is given 

by 

= J±L^ (106) 

However, if we used this vector away from the surface, it would not satisfy ( |l05|) . We could still 



use either expression in ( 104 ) with this vector — we would just have to project the second one 
tangent to the surface. Alternatively, we can use 



(a/-a/)i/2 



(a/.a/)3/2 



(a/-a/)5/2 



f + oif). 



(107) 
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The 0{f'^) terms don't contribute to (104) or (105) on the surface. The contribution of the 0{f) 
terms on the surface to (104) are normal to the surface and ensure that n satisfies (105). 
Either way, on the surface, we get 

(9/-a/)i/2 + (a/-a/)3/2 (5^.^^)5/2 ■ y > 

As this is perpendicular to n, it doesn't matter if we contract its indices with the full metric g^u 
or the induced metric h^,y. We get 

e = e- = BL + ^Wv^ .Qg. 

^ {df-dfy/^ (a/-a/)3/2 ■ ^^"""^ 

B Energy, angular momentum, and entropy of the ring and 
torus in 4+1 dimensions 

The energy, E, angular momentum, L, and entropy S for the 4 + 1 dimensional ring are obtained 



from (|9l|) by differentiation (33) and are given by 



E =- 



(-2379%/6v^^T2i?''/6 - 180V6v/ST2i?i^/^ 



3V6-R + 13(9i? + 26)e 
-1128%/6v/i? + 2i?"/^ + 180V6R + 13i?^ + 1188V6i? + 13i?2 + 2732V6i? + 13i? 

-1690V6v/i?T2^ + 2184V6i?+ 13) 
27r2i?4 



(^-2598156(i? + 2)2/V6i? + 13i?^ 



1354(i? + 2)5/3V6i?+ 13(9i? + 26)3 
-34169688(i? + 2)^/ V6i? + 13i?^ - 186252048(i?+ 2)2/V6i?+ 13i?'' 
-539837568(i? + 2)^/ V6i? + 13i?^ - 879420672(i? + 2)^/ V6i? + UR'^ 
-765123840(i? + 2)^/ V6i? + 13i? 

+A/6(6i?+ 13)(37?(i?(3i?(3i?(243i?(66i?+ 853) + 1105784) + 9364772) 
+44322304) + 37098880) + 38667200) \/R - 278403840(i? + 2)^/3 V6i? + 13^ e 



(468(i? + 2)2/3 V36i?+ 78 



9v/i?T2(9i? + 26)e 

(-396(i? + 2)2/V36i? + 78i?^/^ - 90(i? + 2)^/3 V36i? + 78i?"/^ + 540i?^ 



^2^25/6 ^ 26 

+3276i?2 + 6591i? + 4394)) + ^ ^y,^,^ ^ {uumMOVMRTlSRy^ 

+433026V36i?T78i?3^/'^ + 5596533V36i? + 78R^^/^ + 29856168V36i? + 78i?25/6 
+84282948\/36i? + 78i?^^/^ + 132966912\/36i? + 78i?"/^ - 2598156(i? + 2)^/^7?*^ 
-34169688(i?+ 2)2/^i?,^ - 186252048(i?+ 2)2/3i?4 _ 539837568(i? + 2)2/3i?3 
-879420672(i?+ 2)2/3i?2 _ 765123840(i?+ 2)^/3^ + 38667200V36i? + 78v^ 
278403840(i?+ 2)2/3") e + 0{e^), 
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L = --^ (^TT^elR^VGR + 13 (-6i? + v^fiT2\/36i? + 78v^ - 18 



-1524^36^ + 787?^^/*^ + 1458(i? + 2)^/^R^ + 9630(i? + 2)^/^R'^ 
-21888(i? + 2)2/3i? - 2288V36i?+78v^ + 17160(i? + 2)^/3) + ©(e^), 



These expressions for E, S and L are used for the plot in fig.0. 

For the torus in 4 + 1 dimensions, the energy, entropy and angular momenta can be expressed 
in terms of the derivatives of T In Zgc as in (^J) as well. Using ( |102| ) we find 



E - ^'^o^sin(2xo) / ^1/6(1 ^ i?)i/Vl2 + 9i?(20 + 3i?(8 + 3R)) 

+9(12 + i?(18 + i?(ll + 3i?)))) + !.'^'^^7!)5T3^ (^(^^ + 4)V9i?+T2 (i? 
(9i?(i?(135i? + 554) + 867) + 5564) + 1560) 

-3(i? + lf^^{R{3R{SR{3Ril3bR + 734) + 4907) + 16948) + 30658) + 7904) 
+ 0(e^), 

S = ^'^§-^"^(2X0 ) / ^ lf/^V9R + T2{R{3R + 8) + 8) 
3 + le ^ 

-^(3i? + 4)(3i?(3i?^ + 8) + 14)) - -'R''yR + ^)-^<^Xo) /g^^ ^ ^^./a 

/ 144(V3(i?+ 1)V3) V 

VSR + A{R{SR{R{U5R + 554) + 891) + 2020) + 624) 

- V3v^(3i? + A){R{9R{R{135R + 554) + 861) + 5440) + 1482)) + O (e^) , 



^ 4V2.^^Bi?^nr(Xo)cos^(Xo) /g^^ ^ 2)VwTl ^/3^^i?+T(3i^ + 4) 



3e 

^ .^4fl^y3i? + 4eos(2xo)csc(xo) 3^ ^ 4)^/9i^+12(i?(27i? + 62) + 39) 

36\/2(i? + l)5/3e V ' ' ' 

-3(i? + lf/^{R{3R{27R + 98) + 401) + 208)) + O (e") , 

^ 4V2vr3^3^siV(xo)cos(xo) ^^3^^ ^ 2)^/3i?T4 - V3^^i?TT(3i? + 4)) 
+ (3(i? + l)V3(^(3i?(27i? + 98) + 401) + 208) 



^Ri3R + 4)V9i?+ 12(i?(27i? + 62) + 39)) + O (e°) . 
We use these expressions for E, S, Li and L2 to obtain the plots in fig.[l0|. 

C Notation 

We work with the metric signature ( — !- + + •••)• We use the following index conventions: /x, z/, . . . 
denote space-time indices and a, 6, . . . label planes of rotation, angular velocities, etc. Other 
symbols used are listed below. 
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a'' 
A 

d 
E 
E 

f 

m 
gm 

V 

pfj.u 
Pa 

Qx 



R. 



R 

avg 

r 
S 

s 

S 

T 
f 
T 
% 

t 

Wa 
W 

X 

y{xa) 



acceleration (0) 
surface area of fluid 
m-dimensional ball 
integration constant 
dimensionality of fluid 
total energy of configuration 
see (H) 

surface at / = 
surface at p = g{9) 
terms in e expansion of g (|6: 
induced metric of surface 
proper fluid pressure 
projection tensor (W) 
projector for ( |48| ) 
conserved charge for current X 
heat flux (0) 
radius of fluid tube 
average radius of fluid tube ( [T^ 
radial coordinate (i5) 
radial coordinate (45) 
total entropy of configuration 



mm) 



see ( |36[ ) 

m-dimensional sphere 
proper fluid entropy density 
overall temperature 



see ( |36D 

proper fluid temperature 
transition temperature 
stress tensor 
time coordinate 
velocity vector 
scaled angular velocity (|o|),( 



vector field normal to deformed surface 



see (49) 



see (35) 



surface at r = y{xa) 
grand partition function (| 



tM 

L 
L 
C 

^mfp 
la 

n 

a 

7 
e 

c 

V 



see (36) 



e 

K 

e 
p 

Po 

Pc 

a 
(^s 

X 



-1/2 



surface at r = hira) 
entropy current 
Killing vector for velocity 
total angular momentum 
see (H) 

distance of fluid from origin 

terms in e expansion of C (^,(|62|),(86),(|95|) 

thermalisation scale 

rotational Killing vector 

no. of angular momenta 

unit normal of surface 

equation of state parameter (|]) 

velocity normalisation (1 — v^Y 

expansion parameter (|50|) , (|62|) 

bulk viscosity (^) 

arbitrary Killing vector 

shear viscosity (^ 

angular coordinate (59),( p3|),( 9^), 

or coordinate on (|45) 

expansion (^ 

extrinsic curvature of surface ( |108|) 

trace of extrinsic curvature ( |109D 

thermal conductivity (^ 

surface thickness ( pO| ) 

see (111) 

proper fluid density in §(|2|) 
or shifted radial coordinate ([59|) , (p3|) , (|9 
plasma vacuum energy 
confinement density 
surface tension 
surface energy density 
surface entropy density 
shear tensor (^ 



angular coordinate (45) 
angular coordinate of fluid centre 
terms in e expansion of x 
angular velocity 
see (H) 
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